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REVOLUTION  IN  TURBOMACHINES 
By  John  D.  Stanitz  and  Gaylord  0.  Ellis 


SUMMARY 

A  method  of  analysis  is  developed  for  two- dimensional  flow  on 
general  surfaces  of  revolution  in  turbomachines  with  arbitrary  blade 
shapes.  The  method  of  analysis  is  developed  for  steady,  conqpressible, 
nonviscous,  irrotational  flow  that  is  assumed  uniform  normal  to  the 
surfaces  of  revolution.  Incompressible  solutions  on  a  mean  surface  of 
revolution  between  the  hub  and  shroud  are  presented  for  four  flow  rates 
through  each  of  two  centrifugal  impellers  with  the  same  hub- shroud  con¬ 
tours  but  with  different  blade  spacings.  In  addition,  correlation 
equations  are  developed  whereby  the  velocity  components  and  the  stream 
function  distribution  can  be  predicted  for  compressible  or  incompress¬ 
ible  flow  in  straight-bladed  impellers  only,  with  any  tip  speed,  flow 
rate,  area  variation,  blade  spacing,  and  for  any  flow  surface  of  revo¬ 
lution. 


INTRODUCTION 

In  order  to  achieve  substantial  improvements  in  compressor  and  tur¬ 
bine  efficiency,  new  methods  of  design  must  be  developed,  and  these 
methods  should  be  based  on  detailed  knowledge  of  flow  conditions  within 
turbomachines  (compressors  and  turbines)  .  The  purpose  of  the  analysis 
method,  developed  at  the  NACA  Lewis  laboratory  and  presented  in  this 
report  rs  to  provide  a  tool  whereby  increased  knowledge  of  the  flow 
characteristics  in  turbomachines  can  be  acquired.  In  addition  the  two- 
dimensional  numerical  results  presented  in  this  report  provide  the  basis 
for  future  comparison  with  the  results  of  three-dimensional  solutions. 

Several  two-dimensional  methods  of  analysis  have  already  been 
developed  (references  1  and  2,  for  example).  In  general,  these  various 
methods  differ  in  technique  and  scope  of  application.  In  particular, 
in  reference  2  a  method  of  analysis  is  developed  for  two-dimensional 
compressible  flow  in  turbomachines  with  conic  flow  surfaces  generated 
about  the  axis  of  the  turbomachine  by  straight  lines.  In  this  report 
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the  analysis  method  of  reference  2  is  extended  to  include  arbitrary  sur¬ 
faces  of  revolution  generated  about  axes  of  turbomachines.  In  addition 
correlation  equations  are  developed  for  the  rapid  estimation  of  flow 
conditions  in  impellers  with  straight  blades. 


METHOD  OF  ANALYSIS 
Preliminary  Considerations 

For  given  compressor  operating  conditions,  a  complete  solution  for 
the  flow  through  turbomachines  depends  on  the  machine  geometry  (three- 
dimensional  effects)  and  on  the  fluid  properties  (compressibility  and 
viscosity) .  The  effects  of  compressibility  are  considered  in  the  analysis 
method  of  this  report  but  the  effects  of  viscosity  are  not.  Some  three- 
dimensional  effects  are  indicated  by  the  combination  of  two-dimensional 
solutions.  The  idea  of  these  quasi- three-dimensional  solutions  is  con¬ 
sidered  next. 

Quasi- three-dimensional  solutions.  -  Consider  the  flow  of  an  ideal 
fluid  through  a  typical  passage  between  blades  such  as  shown  in  figure  1. 
The  fluid  is  free  to  follow  whatever  path  the  pressure  and  inertia  forces 
require  of  it.  If,  however,  the  number  of  blades  in  the  turbomachine  is 
considered  to  approach  infinity,  the  space  between  blades,  as  well  as  the 
blade  thickness,  approaches  zero  and  the  path  of  the  fluid  is  limited  to 
the  curved  mean  surface  of  the  blade.  The  fluid  motion  is  thus  reduced 
from  a  general  three-dimensional  motion  to  a  two-dimensional  motion 
restricted  to  the  mean  blade  surface.  The  streamlines  of  this  two- 
dimensional  motion  can  be  projected  on  the  meridional  plane  (axial-radial 
plane)  as  shown  in  figure  2 .  The  solutions  are  axially  symmetrical 
because  the  infinite  number  of  blades  prevent  the  fluid  properties  from 
varying  about  the  axis  of  the  turbomachine. 

Axial- symmetry  solutions  establish  mean  circumferential  flow  con¬ 
ditions,  but  for  finite  blade  spacing  flow  conditions  vary  circumferen¬ 
tially,  as  well  as  from  hub  to  shroud.  In  order  to  investigate  the  cir¬ 
cumferential  variation,  which  is  super inqposed  on  the  circumferential 
mean  flow  conditions  obtained  by  the  axial- symmetry  solution,  it  is 
assumed  that  the  flow  is  limited  to  surfaces  of  revolution  generated  by 
rotating  the  meridional  streamlines  of  the  axial- symmetry  solution 
(fig.  2)  about  the  axis  of  the  turbomachine.  If  these  meridional  stream¬ 
lines  are  selected  vith  stiff iciently  close  spacing,  flow  conditions  can 
be  considered  uniform  normal  to  these  surfaces  of  revolution.  Thus  the 
flow  is  limited  to  two-dimensional  motion  on  flow  surfaces  of  revolution. 

Blade-to-blade  solutions  can  be  obtained  for  every  flow  surface  of 
revolution  generated  by  the  streamlines  in  the  meridional  plane.  There¬ 
fore,  flow  conditions  can  be  determined  approximately  throughout  the 
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passage.  This  resulting  quasi- three- dimensional  solution  is  obtained  by 
the  combination  of  two  types  of  two-dimensional  solutions  axial- symmetry 
solution  and  blade-to-blade  solution.  Although  this  combination  of  two- 
dimensional  solutions  does  not  result  in  a  rigorous  three-dimensional 
solution,  because,  for  example,  the  flow  has  been  constrained  to  surfaces 
of  revolution,  the  combination  must  result  in  a  better  picture  of  the 
true  flow  conditions  than  does  either  of  the  two-dimensional  solutions 
alone.  Analysis  methods  for  the  axial- symmetry  type  of  solution  are 
developed  in  references  3  and  4,  for  example.  In  this  report  analysis 
methods  are  developed  for  the  blade-to-blade  type  of  solution  on  flow 
surfaces  of  revolution  generated  around  the  axis  of  the  turbomachine  by 
known,  or  assumed,  streamlines  in  the  meridional  plane. 

Coordinates .  -  The  cylindrical  coordinates  R,  0,  and  Z  shown 
in  figure  3  are  convenient  to  -use  in  theoretical  studies  of  the  flow  in 
turbomachines.  (All  symbols  are  defined  in  the  appendix.)  These  coor¬ 
dinates  are  dimensionless,  the  linear  coordinates  R  and  Z  having 
been  divided  by  the  blade-tip  radius  r^  (so  that  R  equals  1.0  at  the 

blade  tip,  for  example).  The  coordinate  system  rotates  with  the  blade 
row  about  the  Z-axis.  The  angular  velocity  o>  is  always  considered 
positive  and  in  the  counterclockwise  direction  as  shown  in  figure  3. 


Velocity  components. 


coordinate  system  has  components 
Z  directions,  respectively  (fig.  3) 


The  velocity  Q  relative  to  the  rotating 

Q0,  and  Qz  in  the  R,  0,  and 

These  velocities  are  dimensionless. 


having  been  divided  by  the  stagnation  speed  of  sound  cQ  upstream  of 
the  blade  row,  where 


cn2  =  ygR*T. 


(1) 


In  equation  (1)  R*  is  the  gas  constant,  y  is  the  ratio  of  specific 
heats,  and  T0  is  the  stagnation  (total)  temperature  upstream  of  the 

blade  row.  The  blade-tip  speed  is  likewise  dimensionless  and  becomes  the 
conventional  blade-tip  Mach  number  Mp,  which  is  defined  by 


Mp 


(2) 


Thus  the  tangential  velocity  of  the  blade  at  any  radius  R  is  equal  to 
RMrp  and  the  absolute  tangential  velocity  of  the  fluid  is  equal  to 

(RM^i  +  Qg)  •  -  , 
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In  addition  to  the  three  components  of  the  relative  velocity  Qp, 

Q q,  and  Qg  it  is  convenient  to  define  a  fourth  velocity  component, 

lying  in  the  meridional  plane  (RZ-plane)  and  defined  as  the  meridional 
velocity  Q^,  where  from  figure  3 

2  2  _  2 

Qm  =  +  (3) 

The  meridional  velocity  Qjj  helps  define  the  angle  a,  shown  in  fig¬ 
ure  3,  from  which  figure 

Or  =  s^n  a  (^) 

Qz  =  Qm  cos  a  (5) 

Thermodynamic  relations.  -  From  the  general  energy  equation,  in  the 
absence  of  heat  transfer,  it  can  be  shown  that  along  a  streamline  the 
static  (stream)  temperature  T  is  related  to  the  relative  velocity  Q 

by  (reference  2) 


=  i  +  [(a*4!?)2  *  Q2  -  2Mr*u]  (6) 

where  the  subscript  U  refers  to  conditions  far  upstream  of  the  blade 
row  and  where  A  is  the  whirl  ratio  defined  as  the  absolute  moment  of 
momentum  divided  by  r^cQ  and  given  by 

A  =  R(R%  +  Qe)  (7) 

For  isentrppic  flow  the  pressure  ratio  F  and  the  density  ratio  p/pc 
are  likewise  related  to  the  relative  velocity  Q  by 

X  v  -JC- 

F  =  ^  T  =  ji  +  ^  [(R%)2  -  ft2  -  (e) 

and 

t  ■  (i)""1  'I1  +  ^  [(RMr)2  -  q2  -  2vj) r  1  <9) 

where  p  qnd  p  are  the  static  pressure  and  density,  respectively. 
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Analysis 

A  streamline  that  was  obtained  from  an  axial- symmetry  solution  in 
the  meridional  plane  (fig.  2,  for  example)  is  shown  in  figure  4  together 
with  the  coordinates,  velocities,  and  so  forth  at  a  point  on  the  stream¬ 
line.  At  this  point  a  tangent  to  the  streamline  in  the  meridional  plane 
intersects  the  axis  of  the  turbomachine  with  the  angle  a,  as  indicated. 
If  this  tangent  line  is  rotated  about  the  axis,  it  generates  the  surface 
of  a  cone  that  is  tangent  to  the  flow  surface  of  revolution  obtained  by 
rotating  the  streamline  (fig.  4)  .  At  the  circumf erential  line  of 
tangency  a  fluid  particle  on  the  flow  surface  of  revolution  is  also  on 
the  conic  flow  surface.  This  fluid  particle  is  shown  in  figure  5  on  a 
developed  view  of  the  cone  surface.  The  height  ratio  H  of  the  fluid 
particle  is  defined  by 


where  h  (fig.  2)  is  the  local  spacing  between  center  lines  of  the  chan¬ 
nels  formed  by  the  adjacent  streamlines  (fig.  2)  on  either  side  of  the 
streamline  that  generates  the  surface  of  revolution. 


Absolute  irrotational  motion.  -  In  the  absence  of  shock,  viscous 
effects,  and  so  forth,  the  absolute  motion  of  the  fluid  everywhere  on  the 
surface  of  revolution  is  isentropic  and  hence  irrotational  so  that  the 
absolute  circulation  dF  around  every  fluid  partical  is  zero  and  from 
figure  5 


[(w*p 


Qe)R  ae 


% 

sin  a 


from  which 


[e 


■9 

0  =  2RM]1  sin  a  +  Qg  sin  a  +  R  sin  a  (11) 

Continuity,  *  From  continuity  considerations  for  the  fluid,  particle 
in  figure  5 


(SL  Qg  H  \ 

\Po  sin  a/ 


where  H  is  a  function  of  H  along  the  surface  of  revolution 
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A  stream  function  + 
as 


satisfies  the  continuity  equation  (12)  if  defined 


TUgnfiformation  of  coordinates.  -  For  convenience  in  the  relaxation 
solution  the  following  transformation  of  coordinates  is  introduced: 


ft 


¥ 


d|  =  ■  .  1  . - 55  (14a) 

*  (A 0)  sin  a  R  v 

and 

dT»  -  H  (l4b) 

where  A6  is  the  angular  spacing  of  the  "blades  about  the  axis  of  the  * 

turbomachine.  Equation  (14a)  can  be  integrated  to  give  £  as  a  function 
of  R  with  |  arbitrarily  equal  to  zero  when  R  equals  one,  and  equa- 

Q 

tion  (14b)  can  be  integrated  to  give  r\  equal  to  By  this  trans¬ 

formation,  which  is  conformal,  the  curved  flow  surface  of  revolution 
becomes  a  flat  plane  (Jri-plane).  Lines  of  constant  £  and  T|  are  shown 
transformed  onto  a  flow  surface  of  revolution  in  figure  6.  The  problem 
is  solved  in  the  |i)-plane  and  transformed  bach  to  the  curved  flow  surface 
of  revolution.  The  advantages  of  this  transformation  are:  first,  a 
warping  of  the  blade  shape  -which  for  straight  blades  results  in  straight 
parallel  boundaries  (blades)  in  the  transformed  plane,  thus  permitting 
equal  grid  spacing  for  solution  by  relaxation  methods;  and  second,  a 
reduction  in  the  number  of  terms  in  the  final  differential  equation  for 
the  stream  function  distribution,  thus  reducing  the  labor  involved  in 
obtaining  a  relaxation  solution. 

General  differential  equation  for  stream  function  distribution.  - 
In  terms  of  the  transformed  coordinates,  equation  (11)  for  absolute  irro- 
tational  motion  becomes 


0  =  2RMp  A0  sin  a  +  Qg  A0  sin  a  + 


A 

V 


(15) 
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* 


J. 


* 


I 

* 


and  equation  (13)  for  the  stream  function  becomes 

if 'i 

r  (16) 


Equations  (15)  and  (16)  combine  to  give  (note  that  H  and  R  are  inde¬ 
pendent  of  q) 

2  2  Sin—  Sin  — 

S  ^  S  S  In  H  S\|r  _ Po  d\jr  _ Po  Jhjr 

d£2  +  Sq2  "  “  d|  S|  "  Sq  Sq  _ 


2  -£.  HR2  (A0)2  Mm  sin  a  =  0  (17) 

Po 

P 

In  this  equation  the  density  ratio  —  is  related  to  the  velocity  com- 

Po 

2  2  2 

ponents  and  Q0  by  equation  (9)  where  Q  =  QM  +  Qg  ,  and  the 

velocity  components  are  related  to  the  stream  function  \|r  by  equa¬ 
tion  (16).  This  system  of  equations  ((9),  (16),  and  (17))  can  be  solved 
for  given  boundary  conditions  by  relaxation  methods  (references  2  and  5) 
to  give  the  distribution  of  i|r  for  two-dimensional  flow  past  arbitrary 
blade  shapes  on  the  flow  surface  of  revolution.  The  velocity  conq>onents, 
and  therefore  the  pressure,  density,  and  so  forth,  can  be  obtained  from 
the  distribution  of  i| r  by  means  of  equation  (16)  . 


Incompressible  flow.  -  For  incompressible  flow  the  stagnation  speed 
of  sound  contained  in  the  definition  of  Q  and  is  undefined  and 

its  use  can  be  eliminated  by  considering  the  ratio  which  is  the 

jyirp 

ratio  of  the  relative  velocity  to  the  impeller  tip  speed.  Thus,  with 
—  equal  to  one  for  incompressible  flow  equation  (16)  becomes 


sir =  %  m  (A0) 

% 

sr  =  "  %  m  (Ae) 


(16a) 
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where  from  equations  (16)  is  related  to  +  hy 


*,  =  ± 
yfi  % 


so  that  equation  (17)  becomes  after  dividing  hy  Mrp 
h\  SlnH^i  2  ,aM2  , 

+  '  “5f"  H  "  (Ae)  a  =  9 


(17a) 


which  is  the  partial  differential  equation  used  for  the  numerical  solu¬ 
tions  presented  in  this  report. 


Method  of  Solution 

Numerical  methods.  -  Solutions  of  the  differential  equations  (17) 
or  (l7a)  for  given  blade  shapes  on  given,  arbitrary  surfaces  of  revolu¬ 
tion  and  for  given  flow  rates  and  blade-tip  Mach  numbers  are  obtained  by 
relaxation  methods  using  techniques  discussed  in  reference  2. 

Boundary  considerations.  -  As  in  reference  2  it  can  be  shown  from 
continuity  considerations  that  if  ijr  is  zero  along  the  driving  face  of 
the  blade  (blade  surface  in  the  direction  of  rotation)  the  stream  func¬ 
tion  along  the  trailing  face  is  constant  and  equal  to 

tt  =  'PCAe)  (18) 

In  equation  (18)  the  channel  flow  coefficient  cp  is  defined  by 


PoaTco 


where  W  is  the  flow  rate  through  the  channel  and  aT  is  the  channel 
flow  area  at  the  blade  tip 


(19a) 


For  incompressible  flow 
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The  Joukowski  condition  and  other  boundary  conditions  are  discussed 
in  reference  2. 


MJMERICAL  EXAMPLES 

Solutions  are  presented  for  four  flow  rates  through  each  of  two 
centrifugal  impellers  with  the  same  hub-shroud  contours  but  with  differ¬ 
ent  blade  spacings.  The  solutions  are  for  incompressible  flow  on  a  sur¬ 
face  of  revolution  generated  by  the  mean  streamline  of  an  axial- symmetry 
solution. 


Design  and  Operating  Conditions 

Impeller  geometry.  -  The  hub-shroud  profile  of  the  impellers  is 
described  in  figure  7 .  The  impeller  blades  are  straight  and  thin,  and 
for  mathematical  simplicity  extend  indefinitely  far  upstream  parallel 
with  the  axis  of  thfe  impeller.  The  diffuser  is  vaneless.  Solutions 
were  obtained  for  two  blade  spacings  that,  after  accounting  for  the 
blade  thickness,  measured  14.67°  and  32.80°  about  the  axis. 

Operating  conditions.  -  The  fluid  was  considered  to  be  incompress¬ 
ible  and  nonviscous.  The  solutions  were  obtained  on  a  flow  surface  of 
revolution  generated  by  the  mean  (0.5)  streamline  of  the  axial- symmetry 
solution  in  figure  2  (the  R,Z-coordinates  of  which  streamline  were 
obtained  from  results  of  reference  4  and  are  given  in  table  I  together 
with  sin  a  and  H) .  Four  solutions  for  different  flow  rates  were 
obtained  for  each  of  the  two  blade  spacings.  At  low  flow  rates  an  eddy 
(reference  2)  exists  on  the  driving  face  of  the  blade  and  at  zero  flow 
(closed  throttle  but  impeller  rotating)  the  eddy  occupies  the  entire  pas¬ 
sage.  As  indicated  in  reference  2  the  eddy  flow,  which  consists  of  fluid 
that  remains  in  the  impeller  and  rotates  relative  to  the  impeller  with  an 
angular  velocity  equal  and  opposite  to  that  of  the  impeller,  results 
from  the  condition  of  absolute  irrotational  fluid  motion.  The  four 
solutions  for  each  blade  spacing  were  obtained  for  the  following  flow 
conditions:  (l)  eddy  flow  (no  through  flow  so  that  the  eddy  occupies 

entire  channel) ,  (2)  zero  eddy  flow  (through  flow  rate  just  sufficient 
to  eliminate  eddy) ,  (3)  half  eddy  flow  (half  of  the  through  flow  rate 
required  to  eliminate  eddy) ,  and  (4)  through  flow  rate  used  in  axial- 
symmetry  solution  (fig.  2)  of  reference  4.  These  flow  rates  can  be 
expressed  by  ratios  of  upstream  axial  velocity  to  impeller  tip  speed 

Tke  through  flow  rate  required  to  eliminate  the  eddy  varies 

with  the  blade  spacing  so  that  the  values  of  (Qz/mt)u  for  the  four 

flow  conditions  described  are  given  for  both  blade  spacings  in  the 
following  table: 
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/Qz\ 

ytr/u 

^^-^Blade  spacing 
Flow^^-^A6 y  deg 
condition 

14.67 

32.80 

Eddy  flow 

0 

0 

One  half  eddy  flow 

.060 

.105 

Zero  eddy  flow 

.120 

.210 

Flow  rate  of 

.5429 

.3429 

axial  symmetry 
solution 

The  solutions  for  the  various  ratios  of  flow  rate  to  impeller  tip 
speed  were  obtained  by  superposition  of  the  eddy  flow  solution  (rotating 
impeller  with  no  through  flow)  with  various  percentages  of  the  through 
flow  solution  (flow  through  nonrotating  impeller) . 

A 


Results 

The  results  of  the  numerical  examples  are  presented  in  figures  8 
to  11.  The  results  are  plotted  on  the  surface  of  revolution  between  two 
impeller  blades.  The  |,r|  coordinate  system  is  shown  on  the  surface  of 
revolution.  The  relation  between  |(A0)  and  the  specified  cylindrical 
coordinates  R  and  Z  along  the  surface  of  revolution  was  obtained 
from  a  numerical  integration  of  equation  (14a)  and  is  given  in  table  I. 
From  equation  (14b)  T]  is  related  to  6  by 


Strepml  ines ,  -  The  streamline  patterns  on  the  surface  of  revolution 
are  shown  in  figures  8  and  9  for  the  32,80°  and  the  14.67°  angular  blade 
spacings,  respectively.  For  each  blade  spacing  the  streamlines  are  given 
for  the  four  specified  flow  rates.  The  impeller  rotates  in  the  clockwise 
direction  so  that  the  driving  face  of  a  blade  is  on  the  left  of  the 
channel  and  the  trailing  face  on  the  right.  Except  for  the  complete  eddy 
flow  solutions  (figs.  8(a)  and  9(a))  the  streamlines  are  designated  by  a 
stream  function  ratio  ♦j/OkiJt  suck  that  the  value  of  this  ratio  indi¬ 
cates  the  percentage  of  total  channel  flow  between  the  streamline  and  the 
driving  face  of  the  blade.  (For  the  complete  eddy  flow  solutions  the 
stream  function  ratio  tj/(^j )Tn1n  indicates  the  percentage  of  total 

flow  rate  circulating  within  the  channel  between  the  streamline  and  the 
nearest  blade.)  At  each  radius  the  spacing  of  the  streamline  is  indi¬ 
cative  of  the  relative  velocity,  with  close  spacing  indicating  high 


* 


I 


A 
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velocities  and  wide  spacing  indicating  low  velocities.  Thus,  the  veloc¬ 
ities  are  higher  on  the  trailing  face  of  the  blade  than  on  the  driving 
face,  except  at  the  blade  tip  where  the  velocities  are  equal  on  both 
faces  and  except  for  the  complete  eddy  flow  solutions  (figs.  8(a) 
and  9(a))  where  the  velocities  are  equal  on  both  faces  at  equal  radii. 

As  indicated  in  the  figures  the  fluid  within  the  eddy  remains  in 
the  impeller  and  rotates  in  the  opposite  direction  to  that  of  the  impeller 
so  that  the  velocities  are  inward  along  the  driving  face.  (Also  see 
reference  2,  for  example.)  The  eddy  flow,  which  results  from  the  condi¬ 
tion  of  absolute  irrotational  motion,  is  associated  with  changes  in 
radius  along  the  rotating  flow  surface  and  disappears  as  the  upstream 
surface  becomes  cylindrical  (constant  radius).  This  situation  is  most 
clearly  evident  in  the  solutions  for  no  through  flow  (figs.  8(a) 
and  9(a)).  Thus,  for  quasi- three- dimensional  solutions  the  eddy  flow 
disappears  when  the  surfaces  of  revolution  become  cylindrical.  (This 
condition  is  commonly  assumed  when  two-dimensional  cascade  data  are 
applied  to  the  design  of  axial-flow  turbomachines.)  For  rigorous  three- 
dimensional  solutions  the  eddy  continues  to  exist  but  its  plane  is  normal 
to  the  axis  of  the  turbomachine  and  results  only  in  a  distortion  of  the 
initially  cylindrical  surface.  If  this  distortion  is  neglected  by 
assuming  the  flow  is  constrained  to  surfaces  of  revolution,  the  eddy 
disappears  on  cylindrical  surfaces. 

In  the  solutions  for  no  through  flow  (figs.  8(a)  and  9(a))  the 
greater  magnitude  of  Ol^)min  for  A0  equal  to  32.80°  than  for  A 9 

equal  to  14.67°  indicates  that  the  magnitude  of  the  eddy  flow  rate 
increases  with  increasing  blade  spacing. 


Relative  velocity  components.  -  Lines  of  constant  relative  velocity 
components  are  shown  in  figures  10  and  11  for  the  32.80°  and  14.67°  angu¬ 
lar  blade  spacings,  respectively.  The  meridional  and  relative  tangential 


velocity  components 


(%  %\ 

\M].  %/ 


are  plotted  and  only  the  solutions 


for  the  eddy  flow  are  considered.  (For  solutions  with  through  flow, 


% 

—  remains  the  same  and 

% 


Si 

Mp 

Q 


is  increased  an  amount,  given  in  table  I, 


/  <4Z\ 

that  is  proportional  to  ^  an^  a  function  of  £  (AG)  only.) 


% 

Lines  of  constant  relative  tangential  velocity  ^  are  plotted  in 


figures  10(a)  and  11(a) .  The  values  of 


% 

Mm 


near  the  impeller  tip  are 
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negative  so  that  the  slip  factor  n  (defined  as  the  ratio  of  the  average 
absolute  tangential  velocity  at  the  impeller  tip  to  the  impeller  tip 
speed)  is  less  than  1.0.  The  resulting  slip  factors  for  the  two  blade 
spacings  are 


Blade 

Slip 

spacing 

factor 

A0 

(cLeg) 

14.67 

0.9111 

32.80 

.8142 

Because  —  is  unaffected  by  the  through  flow  the  values  of  a  are 

Mm 

independent  of  tlie  flow  rate  (also  see  reference  6)  .  The  values  of 
—  at  lower  values  of  R  upstream  of  the  impeller  tip  approach  a  max- 

JMm 

imam  positive  value  and  then  decrease  to  zero  as  the  surface  of  revolu¬ 
tion  becomes  cylindrical.  It  is  noted  that  the  maximum  positive  value 

of  —  is  considerably  larger  for  the  wider  blade  spacing. 


Lines  of  constant  meridional  velocity  —  are  plotted  in  fig- 

Mrp 


Q 


ures  10(b)  and  11(b) .  The  values  of  -j—  are  negative  along  the  driving 
face  and  positive  along  the  trailing  face  with  zero  values  along  the 

Si 

center  line  of  the  channel  between  blades.  The  maximum  values  of  77- 

Mji 

occur  on  the  blade  surfaces  near  the  impeller  tip  and  as  the  radius 

Sm 

decreases  —  decreases  and  approaches  zero  as  the  upstream  surface  of 

revolution  becomes  cylindrical.  It  is  noted,  as  for  77-,  that  as  the 

“r 


blade  spacing  increases  the  magnitude  of 


increases 
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CORRELATION  OF  RESULTS 

The  relaxation  solutions  presented  in  this  report  required  consid¬ 
erable  labor  and  it  is  therefore  desirable  to  have  a  more  rapid  method 
of  analysis  to  serve  as  a  useful  engineering  tool.  The  correlation 
equations  to  be  presented  in  this  section  have  been  developed  for  this 
objective.  These  equations  are  limited  to  impellers  with  straight  blades 
on  arbitrary  surfaces  of  revolution,  but  other  less  accurate  methods  of 
analysis  have  been  developed  for  blades  with  arbitrary  shape  (refer¬ 
ence  7,  for  example)  .  The  correlation  equations  to  be  developed  are  an 
extension  of  the  correlation  equations  developed  for  straight  blades  on 
conic  flow  surfaces  of  revolution  in  reference  6.  This  extension  enables 
the  correlation  equations  to  apply  to  arbitrary  surfaces  of  revolution. 

In  this  report,  as  well  as  in  reference  6,  the  correlation  equations  are 
developed  to  predict  the  distribution  of  the  relative  velocity  components 
and  and  the  stream  function  f  throughout  the  impeller  passage. 

The  correlation  equations  are  developed  from  the  fundamental  conditions 
of  continuity  and  absolute  irrotational  motion  and  from  certain  observed 
results  of  the  relaxation  solutions.  These  correlation  equations  predict 
flow  conditions  in  straight-bladed  impellers  with  any  tip  speed,  flow 
rate,  area  variation,  blade  spacing,  and  for  any  flow  surface  of  revolu¬ 
tion.  The  correlation  equations  determine  these  flow  conditions  from 
various  parameters  that  describe  the  design  and  operating  conditions  of 
the  impeller  and  from  tabulated  values  of  any  known  "standard”  solution 
(such  as  that  given  in  reference  6)  . 

In  the  correlation  equations  of  this  report,  and  of  reference  6,  the 
standard  conditions  (obtained  from  any  known  solution,  reference  6,  for 
example)  and  the  nonstandard  conditions  (being  investigated  by  the  corre¬ 
lation  equations)  are  evaluated  at  the  same  values  of  £  and  tj.  The 
effect  of  the  shape  of  the  surface  of  revolution  on  the  fluid  properties 
at  a  given  point  (|,t])  on  the  surface  of  revolution  is  determined  by 
the  effect  of  the  shape  on  sin  a,  which  term  appears  in  the  correlation 
equations,  and  by  the  effect  of  the  shape  on  which  quantity  is 

related  to  sin  a  by  equation  (14a)  . 

Compressible  flow.  -  In  the  previous  development  of  correlation 
equations  the  local  velocity  was  found  to  be,  in  effect,  made  up  of  the 
average  value  of  at  the  given  radius  and  superimposed  velocity  com¬ 

ponents  caused  by  the  rotation  of  the  channel.  The  superimposed  velocity 
components  were  sufficiently  small  not  to  be  greatly  affected  by  compress¬ 
ibility,  which  was  only  considered  to  affect  the  average  value  of  Q^. 

Because  this  effect  on  is  easily  incorporated,  the  correlation  equa¬ 

tions  in  this  report  are  developed  for  compressible  as  well  as  incompress¬ 
ible  flow. 
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Assume  that  the  Qg  component  of  the  relative  velocity  is  corre¬ 
lated  by 


V  =A(Qg)s  (21) 

where  the  superscript  prime  indicates  approximate  values  determined  by 
the  correlation  equations,  where  the  subscript  s  indicates  standard 
values  (given  in  tables  of  reference  6,  for  example) ,  and  where 

(A0)  sin  a 

A  =  [Mp  (A0)  sin  a] “  ^22) 

For  constant  a  equation  (21)  is  identical  to  the  corresponding  corre¬ 
lation  eqution  of  reference  6.  (Note  the  different  definition  of  A6 
in  reference  6.) 

From  equation  (21)  and  using  the  identical  procedure  described  in 
reference  6,  the  correlation  equation  for  becomes 


V  -  <Vav  +  A  [«M  -  («M)avJs  + 

C  (q  )  r  "1 

A(2q  -  1)  i  -|  S  |^(A0)  sin  a  -  (A 0)g  sin  agJ  +  (R  -  Rs)(Mp  A 0)g  sin 

(23) 

■where  R  is  a  known  function  of  £  and  where  from  one-dimensional  con¬ 
tinuity  considerations  (reference  6) 


av 


(24) 


In  equation  (24)  the  average  density  ratio  Pav/pQ  is  obtained  from 

equation  (9)  by  assuming  that  for  straight  blades  the  relative  tangential 
velocity  is  zero  and  the  meridional  velocity  is  equal  to  its  average 
value  so  that 


s  («M>av2 
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Likewise  the  correlation  equation  for  the  stream  function  \|r  becomes 

h  =  11  +  (["Wav  (£  - ")] „ +  c-i  - 1)  i  {-^r1  [m  « - 

(A0)g  sin  agJ  +  (R  -  Rg)  (Mj  A8)g  sin  agj^  (25) 

For  a  complete  eddy  flow  (no  through  flow)  ^  and  (Q^)  av  are  zero  so 

that  equation  (25)  becomes  indeterminate.  If,  however,  both  sides  of 
equation  (25)  are  multiplied  by 


=  —  (Qj  HR  (A0) 

Tt  pQ  v^M'av  v  ' 

which  equality  is  obtained  from  equations  (18)  and  (24),  equation  (25) 
can  be  evaluated  directly  for  \|rT  . 

Incompressible  flow.  -  For  incompressible  flow  the  stagnation  speed 
of  sound  can,  as  before,  be  eliminated  from  consideration  by  expressing 
Q  and  \|r  as  ratios  of  Hp.  Thus  for  incompressible  flow  the  correla¬ 
tion  equations  become 


(21a) 


where 


A*  -  (A8)  sin  a 

[(A0)  sin  a”  ( 


{^m\  RSPavl  .*  1/^ m\  f 

\  Mr r  /  [_  Mj  J+A\V%/[ 


av 


A*  (2n  -  1) 


(i)s  [(A0)  sin  a  -  C*08  sin  aj  + 


(R  -  R_)  (A0) „  sin  a, 


(22a) 


(23a) 
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and 


(R  - 


Representative  values  of 


Rs)(A0) 


sin  a. 


(25a) 


m 

\M^/,  /  * 


and 


♦i 

Wt 


obtained  by 


the  relaxation  solutions  of  this  report  are  compared  in  figures  12  to  14 
with  the  corresponding  values  given  by  the  correlation  equations  [(21a) , 
(23a)  and  (25a)]  using  the  standard  solution  given  in  reference  6. 
Perfect  agreement  corresponds  to  the  45°  straight  line,  and  it  is  seen 
that  the  correlation  equations  are  accurate. 


Effect  of  flow  surface  shape.  -  The  correlation  equation  (23a)  has 
been  used  to  compute  the  velocity  distribution  on  the  blade  surfaces  for 
incompressible  flow  in  impellers  with  straight  blades,  with  constant 
flow  area  along  the  flow  surface,  and  with  various  shapes  of  flow  sur¬ 
faces  of  revolution.  These  flow  surfaces  are  described  by  the  shape  of 
their  generating  lines  (meridional  streamlines)  in  the  RZ-plane.  In 
order  to  demonstrate  the  effect  of  flow  surface  shape,  calculations  have 
been  carried  out  for  generating  lines  that  are  described  by  the  equation 


R 


1 

2 


(2Z  -  l)n  +  1 


(26) 


The  family  of  generating  lines  is  given  by  various  values  for  the 
parameter  n,  and  as  indicated  in. figure  15,  the  surfaces  in  al 1  cases 
start  at  a  radius  R  of  0.5  (impeller  inlet)  and  end  at  a  radius  R  of 
1.0  (impeller  tip) . 

The  resulting  velocity  distributions  on  the  blade  surfaces  are 
plotted  in  figure  16.  As  the  value  of  n  decreases,  the  difference  in 
velocity  on  the  two  blade  surfaces  toward  the  impeller  tip  decreases  nnfl 
in  the  limit  as  n  approaches  zero  the  velocities  become  equal  on  both 
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surfaces .  At  this  limiting  value  of  n  equal  zero  the  flow  is  axial 
(see  generating  line  for  n  equal  to  zero  in  fig.  15)  and  the  veloc¬ 


ities  on  both  blade  surfaces  are  equal  to  the  average  velocity 


V  MT/av' 


which  for  these  examples  is  constant  because  of  the  constant  prescribed 
flow  area. 

As  the  value  of  n  increases  the  difference  in  velocity  on  the  two 
blade  surfaces  toward  the  impeller  tip  increases  and  in  the  limit  as  n 
approaches  infinity  this  difference  approaches  that  for  a  radial  flow 
impeller  (see  generating  line  for  n  equal  to  infinity  in  figure  15). 

If  the  velocity  distribution  curves  in  figure  16  were  extrapolated 
to  lower  values  of  R,  the  curves  would  cross  because,  as  indicated  in 
figure  15,  the  flow  surfaces  that  are  more  nearly  radial  near  the 
impeller  tip  are  more  nearly  axial  near  the  impeller  inlet,  and  vice 
versa. 


SUMMARY  OP  RESULTS 

A  general  method  of  analysis  is  developed  for  steady,  two- 
dimensional,  compressible,  nonviscous,  irrotational  flow  between  arbi¬ 
trary  blade  shapes  on  arbitrary  flow  surfaces  of  revolution  in  turbo¬ 
machines.  Incompressible  solutions  are  presented  for  four  flow  rates 
through  each  of  two  centrifugal  impellers  with  the  same  hub-shroud  con¬ 
tours  but  with  different  blade  spacings.  In  addition,  correlation  equa¬ 
tions  are  developed  whereby  the  velocity  components  and  the  stream  func¬ 
tion  distribution  can  be  predicted  for  compressible  or  incompressible 
flow  in  straight-bladed  impellers  only,  with  any  tip  speed,  flow  rate, 
area  variation,  blade  spacing,  and  for  any  flow  surface  of  revolution. 


Lewis  Flight  Propulsion  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  December  23,  1951 
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AFPEHDIX 

SYMBOLS 


The  following  symbols  are  used  in  this  report: 

A  parameter,  equation  (22) 

A  parameter,  equation  (22a) 

ap  flow  area  at  channel  tip,  equation  (19a) 

co  stagnation  speed  of  sound  upstream  of  turhcmachine,  equa¬ 

tion  (l) 


g 

H 

h 


Mqi 

n 

P 

P 

Q 

% 

QrjQs^Qz 


R,0,Z 


R* 

rT 

T 

W 


gravitational  acceleration 
height  ratio,  equation  (10) 

spacing  between  center  lines  of  channels  formed  by  streamlines 
adjacent  to  streamline  that  generates  surface  of  revolution, 
figure  2 

blade-tip  Mach  number,  equation  (2) 
exponent,  equation  (26) 
pressure  ratio,  equation  (8) 
static  pressure 

relative  velocity  expressed  as  ratio  of  c0 
meridional  component  of  Q,  equation  (3) 

components  of  Q  in  R,  0,  and  Z  directions,  respectively 

cylindrical  coordinates,  dimensionless  (linear  coordinates  R 
and  Z  expressed  as  ratios  of  rT) 

perfect  gas  constant 

blade-tip  radius 

static  temperature 

total  flow  rate  through  channel  (of  height  h)  between  two 
blades 
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a  angle,  figures  3  and  4 

F  dimensionless  fluid  circulation 

X  ratio  of  specific  heats 

T]  transformed  coordinate,  equation  (l4b) 

A 9  angular  blade  spacing  about  axis  of  turbomachine 

X  whirl  ratio,  equation  (7) 

\x  slip  factor 

t  transformed  coordinate,  equation  (l4a) 

p  static  weight  density 

cp  channel  flow  coefficient,  equation  (19) 

\|r  stream  function,  equation  (l3) 

(D  angular  velocity  of  blade  row 

Subscripts: 
av  average 

i  incompressible 

min  minimum 

o  stagnation  condition  upstream  of  turbomachine 

s  standard  solution 

T  blade  tip 

t  trailing  face  of  blade 

U  far  upstream  of  turbomachine 

Superscripts: 

1  estimated  by  correlation  equations 
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TABLE  I  -  CHARACTERISTICS  OF  MEAN  FLOW  SURFACE  OF 
REVOLUTION  USED  IN  NUMERICAL  EXAMPLES 


t  •  A0 

R 

Z 

sin  a 

H 

(Qz/%)u 

(a) 

-1.792 

0.4915 

-0.3394 

0.0082 

3.659 

1.001 

-1.664 

.4920 

-  .2767 

.0140 

3.652 

1.002  , 

-1.536 

.4930 

-  .2135 

.0218 

3.636 

1.004 

-1.408 

.4950 

-  .1502 

.0333 

3.608 

1.008 

-1.280 

.4978 

-  .0870 

.0555 

3.565 

1.015 

-1.152 

.5021 

-  .0232 

.0900 

3.493 

1.027 

;  -1.024 

.5096 

.0411 

.1422 

3.364 

1.050 

-  .896 

.5212 

.1058 

.2222 

3.160 

1.093 

-  .768 

.5408 

.1710 

.3550 

2.858 

1.165 

-  .640 

.5722 

.2348 

.5372 

2.544 

1.237 

-  .512 

.6205 

.2933 

.7312 

2.236 

1.298 

-  .384 

.6890 

.3408 

.8888 

1.901 

1.374 

-  .256 

.7770 

.3737 

.9678 

1.542 

1.502 

-  .128 

.8803 

.3931 

.9913 

1.239 

1.651 

0 

1.0000 

.4052 

.9968 

1.0000 

1.800 

.128 

1.1368 

.4120 

.9996 

.8700 

1.820 

.256 

1.2919 

.4127 

1.0000 

.8664 

1.608 

aValues  of  (Q^/^T^U  ^or  rates  used  in 

numerical  examples  of  this  report  are  given 
in  the  text  under  Operating  conditions. 
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Figure  7.  -  Hub-shroud  dimensions  of  impeller  for  numerical 
examples.  Vaneless  diffuser;  straight  impeller  blades 
extended  indefinitely  far  upstream  parallel  with  axis  of 
impeller. 
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Estimated  meridional -velocity  ratio 
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Figure  13.  -  Comparison  between  relaxation  and  estimated  values  of 
Correlation  equation  (23a). 


Estimated  value  of  stream  function 
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Figure  14.  -  Comparison  between  relaxation  and  estimated  values  of 
Correlation  equation  (25a). 


Radius  ratio 


Trailing  face 
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Figure  16.  -  Velocity- ratio  along  blades  for  various  shapes  of  flow  surface  of  revolution 
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compressible  or  incompressible  flow  in  straight-  compressible  or  incompressible  flow  in  straight  - 

bladed  impellers  only,  with  any  tip  speed,  flow  rate,  bladed  impellers  only,  with  any  tip  speed,  flow  rate, 

area  variation,  blade  spacing,  and  for  any  flow  sur-  area  variation,  blade  spacing,  and  for  any  flow  sur¬ 
face  of  revolution.  face  of  revolution. 
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